CRITICAL BEHAVIOUR IN THE DENSE PLANAR NJL MODEL 
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We present results of a Monte Carlo simulation of a 2+1 dimensional Nambu - Jona-Lasinio model 
including diquark source terms. A diquark condensate (gg) is measured as a function of source 
strength j. In the vacuum phase (gq) vanishes linearly with j as expected, but simulations in a 
region with non-zero baryon density suggest a power-law scaling and hence a critical system for all 
H > fic- There is no diquark condensation signalling superfluidity. Comparisons are drawn with 
known results in two dimensional theories, and with the pseudogap phase in cuprate superconductors. 
We also measure the dispersion relation E{k) for fermionic excitations, and find results consistent 
with a sharp Fermi surface. Any superfluid gap A is constrained to be much less than the constituent 
quark mass scale Eq. 



In the absence of a reliable non-perturbative method 
for calculating the properties of QCD at high densities, 
models of strongly-interacting matter continue to rely 
on results derived from effective theories such as the 
Nambu - Jona-Lasinio (NJL) model Q. Recently such 
approaches have suggested the intriguing possibility of 
superconducting behaviour in quark matter at high den- 
sity as a result of BCS condensation between diquark 
pairs at the Fermi surface |^J^. Even in simplified mod- 
els, however, systematic methods of calculation are hard 
to find, and up to now mostly self-consistent methods 
have been applied Scenarios such as color- flavor 

locking 1^, however, imply that in the high density phase 
the global U(1)b symmetry in the QCD Lagrangian cor- 
responding to conservation of baryon number is sponta- 
neously broken, leading to superfluid behaviour and a 
massless scalar excitation via Goldstone's theorem. A 
description of the fluctuations of these modes and the re- 
sulting long-range interactions may need a more system- 
atic non-perturbative approach; in this Letter we will use 
numerical results from lattice field theory to argue that 
the situation is more complex than hitherto thought. 

The Lagrangian for the NJL model, incorporating a 
baryon chemical potential ^, is 



^-N.iL = + m + fi'yo)i^ - Y - i^lsTtpy 



(1) 



where -0, -0 carry a global isopsin index acted on by Pauli 
matrices t^. The model in 2-f 1 dimensions has been stud- 
ied with /i 7^ using staggered lattice fermions x,X [0- 
the reader is referred to this paper for further details 
of the formulation and symmetries of the lattice model. 
Apart from the obvious numerical advantages of working 
in a reduced dimensionality, the model in this case also 
has an interacting continuum limit for a critical value of 



the couphng 
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making controlled predictions 



in 0, l/g^ = 0.5, rather far from the continuum limit. 
For fi — 0, this results in a ground state with broken chi- 
ral symmetry and a dynamically generated fermion mass 
So — 0.71, together with a triplet of Goldstone pious. 
The "constituent quark mass" Sq defines the physical 
scale in cutoff units. As fi is raised from zero, the system 
remains essentially unaltered until /i = Hc — 0.65, where- 
upon the constituent quark mass S falls to zero as chiral 
symmetry is restored in a strong first-order transition, 
and baryon density n — {'ijj^o'4') jumps discontinuously 
from zero to a non-zero value; by /i = 0.8, the density 
n ~ 0.25 quarks of each isopsin component per lattice 
site. 

The tendency for diquark pairing at high density was 
investigated via measurements of the diquark propaga- 
tor G{t) = ^^{qq{0,O)qq{t,-x.)) [Q. In particular it was 
found that for fj, > fic and the choice of isoscalar diquark 
operator qq — x*^^ {x)t2x{x) (the "spectral scalar" chan- 
nel of 0), G{t) plateaus at a non-zero value as t — > oo, 
implying {qq){qq) 7^ by the cluster property. Stud- 
ies on varying spatial volumes, however, reveal that the 
plateau height is not an extensive quantity as naively 
expected, and that diquark condensation signalling a su- 
perfluid state is not unambiguously indicated. 

The situation can be clarified by the introduction of 
source terms for diquark and anti-diquark pairs: 

C[jj] = CNjL+Jx'''{x)T2x{x)+Jxi^)r2f'-ix). (2) 

The functional integral may now be written using the 
Gor'icov representation M as 



Z[j,j] = {Pi{A[j,j])), 



(3) 



where (. . .) denotes averaging with respect to bosonic 
auxiliary fields a and 7? introduced to make Cnjl bilinear 
in X aud x, and the antisymmetric matrix A is 



possible in principle. In the simulations of this Letter, as 



A = 



JT2 



jT2 



(4) 
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where M = M(/i; cr, 7?) is the conventional kinetic opera- 
tor oi Cnjl 0- For convenience we now define diquark 
operators qq± via 

qq±{x) = \ W{^)r2X{x) ± x{^)r2t'' i^)] , (5) 

with corresponding sources i± = j ± j. The diquark 
condensate is readily calculated to be 



, , 1 c»lnZ 1 , , 



X 



Xe 



(6) 



Under U(l)s transformations x 
the operators qq± rotate into each other; condensation of 
{qq+) ^ implies a massless Goldstone mode in the qq^ 
channel as 0, illustrated by a Ward identity for the 

susceptibility x-'- 



X- 



J-=0 



{m+) 
i+ 



(7) 



Numerical simulation using the functional measure (j^) 
requires us to address the issue of calculating the Pfaf- 
fian We begin by observing that Pf^ = idet^^, 

the sign depending on the details of the ordering of the 
Grassmann integration. Now, using the property of a 
square block matrix 

det ^ ^ Z ) " detXdet(Z - WX~^Y), (8) 

together with the property T2MT2 = M* , we find det^ = 
det(jj + M'l'M/4), and hence is real and positive if jj is 
chosen real and positive. It follows that Pf^ is real. To 
avoid the sign problem, however, we choose to simulate 

with measure (iei^{A^A) = Pf^^, corresponding to a 
doubling in the number of fcrmion species. This is usual 
practice in simulations of four-fcrmi models; "crosstalk" 
between the two species, which might conceivably lead 
to problems for /i ^ due to light unphysical baryonic 
states in the spectrum, is subleading in this case and 
causes no problems p^ . It is straightforward to imple- 
ment the simulation using a hybrid molecular dynamics 
'R' algorithm — the resulting dynamics are those of 
2 flavors of staggered lattice fermion, corresponding to 
Nf = 4 flavors of four-component physical fermion. 

We have performed simulations on lattice sizes 16"^, 
16^ X 24, 24^, 32^ and 48^, with chemical potentials 
/i = 0,0.2 corresponding to the low density chirally- 
broken phase, and ^ = 0.8 in the high density chirally- 
restored phase Q. Diquark source values ranging from 
j = 0.3 down to as low as j — 0.025 were used; through- 
out we took ] — j and kept a bare Dirac mass m = 0.01 
to assist identification of chirally broken and restored 
phases. We plot ln(gq+) vs. In j in Fig. |l|. At low density 
the results for /i = and /i = 0.2 are very similar and sup- 
port a linear relation {qq+) oc j, suggesting that U(1)b is 
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FIG. 1. In(gq+) vs. Inj 

unbroken in the zero-source limit. For ^ — 0.8, when chi- 
ral symmetry is restored, the situation is different: whilst 
there is still no evidence for a non-zero condensate as 
j ^ 0, the log-log plot now displays a marked curva- 
ture, together with evidence for significant finite volume 
effects. Empirically we find that on a x Lt system 
the dominant correction scales as 1 /Lt - this enables an 
extrapolation of the data with Lt > 24 to the thermody- 
namic limit, shown in Fig. ^. 
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FIG. 2. Extrapolation of {qq+) at /j. = 0.8 to the infinite 
volume limit. 

The extrapolated data, denoted by filled squares in 
Fig. 1^, suggest the best description in the high-density 
phase is a power law: 

{qq+)^r, (9) 

with a ~ 0.31. There is thus no superfluid condensate 
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in this region of parameter space, and no consequent 
Goldstone pole. The Ward identity (0), however, im- 
plies that the transverse susceptibility x- diverges as 
thus long-wavelength fluctuations remain impor- 
tant, and there are strongly-interacting massless states 
present. This is the first main result of this Letter. 

Now, the value /i — 0.8 was originally chosen merely 
as characteristic of the high-density phase; moreover data 
from partially quenched simulations ||^ suggest that the 
power-law scaling @ may be generic for ji > ^c- This 
is analogous to the behaviour of the 2d X-Y model [ISj, 
which in accordance with well-known theorems does 
not have spontaneous magnetisation, but instead exhibits 
a line of critical points for all T < Tc In this low- 

temperature phase the modulus of the spin variables is 
constant, but the phase fluctuates; the divergence of X- 
signals correlation of the phase over macroscopic distance 
scales without long-range phase coherence. This critical 
behaviour without an order parameter can be charac- 
terised by two exponents 6 and rj, defined in our notation 
by {qq+) (xji, and \imj_^^n{qq~{0)qq-ix)) oc l/x'^^'^+'K 
For the X-Y model ^ > 15 and 77 < i (the in- 

equalites being saturated up to logarithmic corrections 
at T = Tc), and satisfy the hyperscaling relation 



it propagates as C{x) (x x' 



, implying massive 



d + 2-r] 
d-2- 



(10) 



For the dense planar NJL model, by contrast, our results 
suggest 5 ~ 3. 

The phenomenon of long-range phase coherence being 
wiped out by soft transverse fiuctuations is particularly 
interesting for a composite order parameter field, sug- 
gesting as it does that some dynamical pairing mecha- 
nism is in play even in the absence of symmetry breaking. 
It has been suggested that such a state may describe a 
pseudogap phase, characterised by a suppression of spec- 
tral weight in the vicinity of the Fermi surface, and ob- 
served as a precursor to the superconducting state for 
underdoped cuprate superconductors It is intrigu- 

ing in the current context that cuprate superconductivity 
is a planar phenomenon |^. 

A field theoretic example of possible relevance is 
studied in |ll7| , where results from an exactly soluble 
2d fermionic model are generalised to the Gross-Neveu 
model with U(l) chiral symmetry [|l^. Once again 
due to the low dimensionality no long-range ordering 
via (ipip) 7^ is possible; however for sufficiently large 
number of fiavors Nf the model has a phase analogous 
to the low temperature phase of the X-Y model, with 
T] ~ 0{1/Nf). The physical fermion is a superposition 
of positive and negative chirality states and has zero net 
chirality; despite the absence of chiral symmetry breaking 



propagation as x — > 00, but also a non-trivial spectral 
function, signalled by the departure from canonical scal- 
ing of the pre-exponential factor. 
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*An important distinction is that in cuprates the pairing 
operator is d-wave, whereas for the NJL model it is s-wave. 



FIG. 3. Dispersion relation E{k), for gapped (left) and 
ungapped (right) cases. 

These considerations motivate a study of the spectrum 
of spin-i excitations in the NJL model by analysing the 
fermion propagator C(k, t) — Re(x(0, 0)x(i, x)e*'' ''), 
where for staggered fermions the sum only includes x 
which are an even number of lattice spacings from the 
origin in each direction. Note that for j ^ 0, the most 
general definition of C{t) contains both "normal" (xx) 
and "anomalous" (xx) components. In this Letter we 
analyse only the isosinglet normal component, with k 
oriented along a lattice axis. Fig. || shows possible spec- 
tra E{k) for /i > /ic, n > 1^. Solid lines denote states 
which are occupied in the ground state; following com- 
mon usage we label these "Dirac" and "Fermi" seas (note 
there is no quantum number which distinguishes between 
these branches). Dashed lines denote empty states, and 
hence describe the spectrum of allowed excitations. Note 
that for k < kp the lowest excitations vacate states in 
the Fermi sea, and hence are 'hole-like': in contrast for 
k > kp, excitations add quarks to the system, and are 
'particle-like'. On the right hand side the curve is contin- 
uous as it crosses the axis; there is a sharp Fermi surface 
with well-defined Fermi momentum k — kp, where exci- 
tations cost zero energy. On the left hand side we have 
sketched the effect of a discontinuity at the Fermi sur- 
face arising from a BCS instability; the lowest excitation 
energy now corresponds to the gap A. 

We have measured C(k, <) on 32"^ lattices at /i = 0.8. 
Empirically we find that for t even it is consistent with 
zero, and that for t odd it is well described by a 3- 
parameter fit C(k, i) = Ae"^* + Be"^(^*"*l For smaU 
values of k we find A ^ _B, ie. the signal predominantly 
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FIG. 4. E{k) at fi = 0.8. 

forwards-moving, and vice-versa for k as approaches its 
maximum value tt/2. As j is decreased the switch be- 
tween these two behaviours becomes increasingly pro- 
nounced. On comparison with the analytic form for free 
fermions, we label these cases as 'hole' and 'particle' re- 
spectively. We also find that E{k,j) decreases monoton- 
ically and roughly linearly with j, and a linear extrapo- 
lation to j = can be made based on j G [0.1,0.025]. 
The results, shown in Fig. ^, are in good qualitative 
agreement with the expectations of Fig. ^ Conserva- 
tively, we can bound any possible gap from above, and 
find A < 0.05 <C So 0.7, in contrast with expec- 
tations based on the self-consistent gap equation ap- 
proach of [||. More boldly, we plot a curve of the form 
E{k) = \P + Q{^ — k)'^\, which suggests it is plausible the 
dispersion relation is continuous near k = kp, and hence 
A = 0. The fit yields a kp — 0.90 slightly larger than the 
Fermi energy /i, which is characteristic of a system with 
attractive inter-particle forces. We conclude that there is 
a sharp Fermi surface in this region of parameter space: 
this is the second main result of the Letter. 

In the future, it will be helpful to explore other values 
of fi > /ic, to see whether d — 6{fi). Independent esti- 
mates of rj from simulations at j = would also be valu- 
able; one could then use hyperscaling (|o|) to probe the 
effective dimensionality of the system. It has been sug- 
gested that absence of diquark condensation in the planar 
NJL model is an artifact due to the low dimensionality 
of the Fermi surface in 2-1-1 dimensions 0; it should be 
mentioned, however, that partially quenched data taken 
on small 3+1 dimensional systems also appear to scale 
according to (||) ||]. Finally, by analogy with cuprate 
superconductivity, it may be the case that this critical 
phase we have found is a precursor to a superfluid phase 
setting in at higher densities. Simulations with the cur- 
rent parameters become plagued by lattice artifacts for 
/i ^ 1.0 ||ll[]. To observe superfluid behaviour simulations 
closer to the continuum limit may be necessary. 
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